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We study a particular class of moving average processes which possess a property 



^ , Abstract 

> 

\^ \ called localisability. This means that, at any given point, they admit a "tangent 

r^ ■ process" , in a suitable sense. We give general conditions on the kernel g defining the 

moving average which ensures that the process is localisable and we characterize the 
nature of the associated tangent processes. Examples include the reverse Ornstein- 



QQ ' Uhlenbeck process and the multistable reverse Ornstein-Uhlenbeck process. In the 

^D ■ latter case, the tangent process is, at each time t, a Levy stable motion with stability 

index possibly varying with t. We also consider the problem of path synthesis, for 
which we give both theoretical results and numerical simulations. 



1 Introduction and background 

In this work, we study moving average processes which are localisable. Loosely speaking, 
this means that they have a well-defined local form: at each point, they are "tangent" to 
a given stochastic process. 

Localisable processes are useful both in theory and in practical applications. Indeed, 
they provide an easy way to control important local properties such as the local Holder 
regularity or the jump intensity. In the first case, one speaks of multifractional processes, 
and in the second one, of multistable processes. Such processes provide fine models for real 
world phenomena including natural terrains, TCP traffic, financial data, EEG or highly 
textured images. 

Formally, a process Y{t) defined on M (or a subinterval of M) is h- localisable at u if 
limr^o+ r~^(Y{u + rt) — Y{u)) exists as a non-trivial process in t for some h > 0, where 



the convergence is in finite dimensional distributions, see [31 H]. When convergence occurs 
in distribution with respect to the appropriate metric on C(M) (the space of continuous 
functions on M) or on D{M.) (the space of cddlag functions on M, that is functions which 
are continuous on the right and have left limits at all t G M), we say that Y is strongly 
localisahle. The limit, denoted by F^ = {Ylit) : t G M}, is called the local form or tangent 
process of y at m and will in general vary with u. A closely related notion is that of locally 
asymptotically self-similar processes, which are for instance described in [2j. 

The simplest localisable processes are self-similar processes with stationary increments 
(sssi processes); it is not hard to show that an sssi process Y is localisable at all u with 
local form FJ = Y . Furthermore, an sssi process Y is strongly localisable if it has a version 
inC(M) or D(M). 

In |5] , processes with prescribed local form are constructed by "gluing together" known 
localisable processes in the following way: let U be an interval with u an interior point. 
Let {X(t,f) : (t,t>) G U x U} be a random field and let Y be the diagonal process 
Y = {X{t, t) : t G U}. In order for Y and X[-, u) to have the same local forms at m, that 
is y^(-) = X'^{-,u) where X^(-,-u) is the local form of X(-,m) at m, we require 

X(M + rt,M + rt) -X(m,m) fdd „,., . ,^ ^. 
-^ ^A„(t,uj (l.lj 

as r \ 0. 

This approach allows easy construction of localisable processes from "elementary 
pieces" which are known to be themselves localisable. In particular, it applies in a straight- 
forward way to processes X(t, v) such that X(-, v) is sssi for each v. 

In this work we shall study a rather different way of obtaining localisable processes. 
Instead of basing our constructions on existing sssi processes we will consider moving 
average processes which, as we shall see, provide a new class of localisable a-stable pro- 
cesses. 

The remainder of this paper is organized as follows: in section [21 we give general condi- 
tions on the kernel g defining the moving average process to ensure (strong) localisability. 
Section [3] specializes these conditions to cases where explicit forms for the tangent process 
may be given, and presents some examples. In section HI we deal with multistable moving 
average processes, which generalize moving average stable processes by letting the stabil- 
ity index vary over time. Finally, section [5] considers numerical aspects: for applications, 
it is desirable to synthesize paths of these processes. Using the approach developed in 
|12j . we first explain how to build traces of arbitrary moving average stable processes. 
In the case where the processes are localisable, we then give error bounds between the 
numerical and theoretical paths. Under mild additional assumptions, an 'optimal' choice 
of the parameters defining the synthesis method is derived. Finally, traces obtained from 
numerical experiments are displayed. 

2 Localisability of stable moving average processes 

Recall that a process {X{t) : t G T}, where T is a subinterval of M, is called a-stable 
(0 < a < 2) if all its finite-dimensional distributions are a-stable, see the encyclopaedic 
work on stable processes pi]]. 2-stable processes are just Gaussian processes. 

Many stable processes admit a stochastic integral representation. Write S'a(o", /?, /x) 
for the a-stable distribution with scale parameter a, skewness /? and shift-parameter /i; 



we will assume throughout that (3 = fi = 0. Let {E,S,m) be a sigma-finite measure 
space {m will be Lebesgue measure in our examples). Taking m as the control measure, 
this defines an a-stable random measure M on ii^ such that for A E S we have that 
M{A) ~ S'a(m(y4)^/", 0, 0) (since /? = 0, the process is symmetric). 
Let 

JFq, = J^a{E,£,m) = {/ : / is measurable and ||/||a < C)o}, 

where || \\a is the quasinorm (or norm if 1 < a < 2) given by 

\ 1/" 
\f{x)rm{dx)] . (2.1) 

E / 

The stochastic integral of / G J-'a{E, S, m) with respect to M then exists [TOl Chapter 3] 
with 

[ f{x)M{dx) ^ S^iaf, 0,0), (2.2) 

Je 

where a/ = ||/||a. 

We will be concerned with a special kind of stable processes that are stationary and 
may be expressed as moving average stochastic integrals in the following way: 

Y{t)= I g{t-x)M{dx) (tGM), (2.3) 

where g G J^a is sometimes called the kernel of Y . 

Such processes are considered in several areas (e.g. linear time-invariant systems) and 
it is of interest to know under what conditions they are localisable. A sufficient condition 
is provided by the following proposition. 

Proposition 2.1 Let < a < 2 and let M he a symmetric a-stable measure on M with 
control measure Lebesgue measure C Let g E J^a and let Y be the moving average process 



Y{t) = g{t- x)M{dx) {t G R). 

Suppose that there exist jointly measurable functions h{t, .) G J-'a such that 

g{r{t + z)) - g{rz) 



lim 

r-»0+ 



7-7 



h{t,z) 



dz = (2.4) 



for all t E M, where 7 + 1/a > 0. Then Y is {'y + 1/ a) -localisable with local form 
rj = {J hit, z)M{dz) : t G M} at allue M. 

Proof. Using stationarity followed by a change of variable z = —x/r and the self-similarity 
of M, 

Y{u + rt)-Y{u) = Y{rt)-Y{0) 

{g{rt — x) — g{—x))M{dx) 



r^/" hg{r{t + z))-g{rz))M{dz), 



where equalities are in finite dimensional distributions. Thus 

Y(u + rt)-Y(u) /",, ,,,,,, f fg(r(t + z))-g(rz) ,, ,\,.,r^ 
^ ^Jyo. - J h(t, z)M{dz) = j [^-^ ir^-^ - ^^^' ^V ^(^^)- 

By [ini proposition 3.5.1] and (El]), r-''-^/''{Y{u + rt) - Y{u)) -^ Jh{t,z)M{dz) in 
probability and thus in finite dimensional distributions. ■ 

A particular instance of (12.31) is the reverse Ornstein-Uhlenbeck process, see [TIH Sec- 
tion 3.6]. This process provides a straightforward application of proposition 12. 1[ 

Proposition 2.2 (Reverse Ornstein-Uhlenbeck process) Let A > and 1 < a <2 and let 

M be an a-stable measure on M with control measure C The stationary process 



/oo 
exp(-A(a; - t))M{dx) {t E 



has a version in D(M.) that is 1/a-localisable at all u E R with Y^ = La, where Lait) := 
Jp M{dz) is a-stable Levy motion. 

Proof. The process y is a moving average process that may be written in the form (12. 3 p 
with g{x) = exp(Ax)l(_oo,o](a^)- It is easily verified using the dominated convergence 
theorem that g satisfies ( 12. 4p with 7 = and h(t, z) = — l(_t,o](^) for t > and h{t, z) = 
— l(o,_t](-2) for t < 0, so proposition 12.11 gives the conclusion with Y^{t) = — M([— 1,0]) = 
La{t) for t > and a similar formula for t < 0. ■ 

Proposition 12.11 gives a condition on the kernel ensuring localisability. With an addi- 
tional constraint we can get strong localisability. First we need the following proposition 
on continuity. 

Proposition 2.3 Let < a < 2, g E J-'a and let M be an a-stable symmetric random 
measure on M with control measure C. Consider the moving average process defined by 
( (^.g|) . Suppose that g satisfies, for all sufficiently small h. 



\g{h — x) — g{—x)\'^dx < c\h\ 



where c > and A > 1. Then Y has a continuous version which satisfies a 6 -Holder 
condition for all 9 < {X ~ l)/a. 

Proof. By stationarity, 

Y{t)-Y{t') = Y{t-t')-Y{0) 

= I {g{t-t' -x)-g{~x))M{dx). 

So for < p < a 



E|F(t)-r(t')r < cJj\g{t-t'-x)-g{-xTdx] 



p/a 



< C2|t-t'|^P/". 

The result then follows from the Kolmogorov criterion by taking p arbitrarily close to a. 



Proposition 2.4 With the same notation and assumptions as in proposition \2.1\ suppose 
that in addition that g satisfies, for all sufficiently small h, 

f \g{h -x)- ^(-x)r dx < c\hr~^\ (2.5) 

where c > and'y > 0. Then Y has a version in C{M.) that is {'j + l/ a) -strongly localisable 
with Y^ = {/ h{t, z)M{dz) : t G M} at all u e M. 

Proof. By proposition [231 ^(^) has a continuous version and so Zy.{t) := r^('''+^/"^(y(rt) — 
Y{fS)) also has a continuous version. Thus, for < p < a, by stationarity and setting 
h = r\t — t'\ sufficiently small, 



E\Zr{t) - Zr{t')\P = E\Zr{t - t') - ZriO)\ 



V 



Cl 



Cl 



g^r{t-t')-x)-g{-x)'' ^^/" 



g{h - x) - g{-x)\- V^" 



dx 



dx\ |t_t'|(7«+lW" 



provided |t — t'| is sufficiently small, using (12.51) in the last step. We may choose p suf- 
ficiently close to a so that (70; + l)p/« > 1. By a Corollary to Kolmogorov's criterion 
(see e.g. [9l Theorem 85.5]) the measures on C{K) underlying the processes Zr are con- 
ditionally compact. Thus convergence in finite dimensional distributions of Z^- to y^ as 
r \ implies the convergence in distribution (with F^ necessarily having a continuous 
version). Together with localisability which follows from proposition 12. II this gives strong 
localisability. ■ 

Note that the reverse Ornstein-Uhlenbeck process is a stationary Markov process which 
has a version in D{M) see |1T1 Remark 17.3]. It also satisfies (12.51) for a > 1 with 7 = 0. 
However, we cannot deduce that it is strongly localisable since proposition 12.41 is only 
valid for 7 > 0. The case 7 = would be interesting to deal with, but is much harder 
and would require different techniques. 

3 Sufficient conditions for localisability and examples 

For the reverse Ornstein-Uhlenbeck process, it was straightforward to check the conditions 
of proposition 12. 1[ In general, however, it is not easy to guess which kind of functions g 
in Ta will satisfy (12.41) . In this section we will find simple practical conditions ensuring 
this. 

Recall that the following process is called linear fractional a-stable motion: 

/oo 
fa,H{b^,b',t,x)M{dx), 

where t G M, &+, &" G M, and 

/„,h(&+, b-,t, x) = 6+ ((t - x);-^/° - (-x);-^/°) 

+ r((t-x)^-'/"-(-x)^-'/"), (3.1) 



where M is a symmetric a-stable random measure (0 < a < 2) with control measure 
Lebesgue measure. Being sssi, La,H,b+,b- is locahsable. In addition, it is strongly localis- 
able when H > 1/a, since its paths then belong to C(M). 
Recall also that the process 

L^(t) = f M{dx) (3.2) 

Jo 

is a-stable Levy motion and the process 

r+oo 

Za{t)= {ln\t-x\-ln\x\)M{dx) (3.3) 



is called log-fractional stable motion. 

We are now ready to describe easy-to-check conditions that ensure that propositions 
O and [231 apply. 

Proposition 3.1 Let < a < 2, g ^ T^ o,nd M he an a-stable symmetric random 
measure on M with control measure C. Let Y be the moving average process 

Y{t)= I g{t-x)M{dx) (teM). 

// there exist c^, C(^, 7, a, c, ?7 G M with c > 0, r^ > and < 7 + l/a < a < 1 such that 

9{r) + , g{-r) 
> Cn and > On 

as r \0 and 

\g{u + h) - g{u)\ < clhl^lu]^-" {u e R, \h\ < ri) , (3.4) 

then y zs (7 + l/a)-localisable at all u &M. with local form 

(6) FJ = (c+ - Co-)L, z/7 = 0. 

//, in addition, 7 > and < a < 2 then Y has a version in C(M) and is strongly 
localisable. 

Note that condition (13.41) on the increments of g may be interpreted as a 2-microlocal 
condition, namely that g belongs to the global 2-microlocal space Cj'""^, see [7j. Remark 
also that, in order for this condition to be satisfied by non-trivial functions g, one needs 
a < 1, which in turns implies that 7 < 1 — 1/a and — 76 (1/a, 1 — 7]. 

Proof (a) We have 

g{r{t + z))-g{rz) ^(^1^ + ^1) , , , ,71 , 9{-r\t + z\) 

- \t + z\ 'l{t+z>o} + / 1, I — 7r:—\t + z\'l{t+z<o} 



r"/ {r\t + z\p' ' ''^-' {r\t + z\p 

9{-r\z\) 

-{z>0} / I |x^ FI '1{2<0}- 



9ir\z\) gi-r\z 

( I \\A^ -'-{2>0} / I IN 



As r ^ 0, 

g{r{t + z)) - gjrz) 



cl\t + 2;pl{j+^>o} + Cg |t + ^pl{i+2<o} - c,[|2;|^l{^>o} - Cq |zpl{^<o} 
c+(t + z)l - 4(^)T + Co-(t + ^)I - Co^(.)I (3.5) 



To get convergence in L" we use the dominated convergence theorem. Fix e > and 
m > such that for all < m < e, 



9\u) 



u-y 



< m and 



9K.-U) 



u-y 



-Cn 



< m. 



For fixed t write fr{z) = r '^{g{r{t + z)) — (^(rz)). If t = 0, 

Jr\Z) = Ja,'y+l/a{CQ ,Cq ,t, Z) = U, 

thus fr — > fa,'y+i/a{cQ,CQ,t,—.) which belougs to L°'. Assume now t G M*. There is a 
constant nii such that |/j,(z)|" < mi(l + l^l"^ + |t + z|''')° for all \r\ < j-rg and \z\ < 1. 
From (I33D 

for |r| < ri/\t\, so, as (7 — a)a < —1 and 7a > —1, 

mi(l + |zp + \t + z[<Ydz + /" |2|(^-'^)"dz < 00. 

|2|<1 •J\z\>^ 

Since also /a,7+i/a(C(J', c^j", t, — 2;) G L", the dominated convergence theorem implies 
that fr{z) — > /a,7+i/a(c^, 0(5", t, — z) in L". The conclusion in case (a) follows from propo- 
sition [2]T], (13.11) . and noting that M is a symmetric a-stable measure. 

(b) In this case the limit (13. 5p is 



g{r{t + z)) - g{rz) 



7-7 



(c+-Co-)l[o,t](-^) if t>0 
(c+-Co)l[t,o](-^) if t<0 



Dominated convergence follows in the same way as in case (a) so the conclusion follows 
from proposition 12.11 and (13.21) . 

Moving to strong localisability, for h small enough, 

I \g{h-x)-g{-x)\''dx<ci f |/i|"Trfx < C2|/ir^+\ 

J\x\<3\h\ J\x\<3\h\ 



and 



\gih-x)-gi-x)\''dx < ci\hr Ix]^^-"^" dx 

\x\>3\h\ J3\h\ 



and the conclusion follows from propositions 12.31 and [2 

We now give an alternative condition for localisability in terms of Fourier transforms. 
Note that the Fourier transform fa,H{b~^, b~, t, ^) of fa,H{b^, b~ ,t, .) is given by 

lMb^,b-,t,0 =TiH+l - l/a)^^A;i- 

(ITT \ / ITT 

ysgn(0(i^ + 1 - 1/a) j + r exp (^ - —sgn{0{H +1-1/ a 

Proposition 3.2 Let 1 < a < 2, andY be defined by Ii2.3\) . If there exist I = li+il2 E C* , 
7 e (-^, 1 - ^), a G (0, 1 - (7 + ^)) and K G UiR) with p G [1, 1/(7 + ^ + a)), such 
that for almost all ^ > 0, 

e+%) = / + ^%), (3.6) 

then Y is {'-f + l/a)-localisable at a// u G M with local form 

(a) y;; = iv„^^+i/a,6+,6- ifi^o 

{b) Y:^ = H,Z^ + hL^ ^/7 = 0, 
where 



2r(7 + 1) Vcos(7r(7 + l)/2) sin(7r(7 + l)/2) 



2r(7 + 1) Vcos(7r(7 + l)/2) sin(7r(7 + l)/2) 
Proof (a) First note that, with 6"'' and b~ as above, we have, for 2; 7^ 

Set /r(-2) = r^^{g{r{t + z)) — g{rz)). Then /^ G JF^, and 

».7+l "\j,' 



fr(0 = h^9Q- 



now 



With a' such that — I — ^ = 1 we have fr G J-'a' and fa,-y+i/a{b'^,b ,t,^) G JFq,/. We 
show that ll/r — fa,'y+i/a{b~^,b~ ,t,—-)\\a -^ whcu r — >• 0. Note that (13. 6p implies that 
for ^ < 

Writing /(^) = fa,-y+i/a{b'^, b~,t, — ^), for almost all ^ G K 

(e'«' - 1) r"_^ { 






r 



lei 



7+1+ 






Let Hr{u) = K{ru). Then Hr{^) = -K{^) and we may write for a + 7 7^ 

^(0-/(0 = r"+i/,,,+i/„+„(6,6,t,-0^r(0, (3.7) 

where b = l/(2r(7 + a + 1) cos(7r(7 + a + l)/2)). 

It is easy to verify that fa,'y+a+i/a{b'^ , b~ ,t, — ) G L^ for all /3 > 1/(1 — 7 — a). By the 
conditions on a and p, there exists such a 3 which also satisfies - + 1 = 4 + - and in 
particular, - + 4 > 1- Consequently we may take the inverse Fourier transform of (13. 7p 
see, for example, [131 Theorem 78] to get: 



friz) - f{z) = r«+Va,7+l/a+a(&+,&",t, "O * ^r(^) 

where * denotes convolution. As ^ + 1 = 4 + -, the Hausdorff- Young inequality yields 



Wfr — fcx,-y+l/a{b ,b ,t, — )||q < r" || /Q,7+l/a+a ||/3||-f^) 



r lip 



^ a+1 II r II II rxll 

We conclude that fr -^ fa,--/+i/a{b'^,b~' ,t, — ) in L". The result follows from proposition 
12. 1[ The case a + 7 = is dealt with in a similar way. 

(b) Let zt and k be defined by 

'^ ^ \ -l]t,o[{x) ift<0 

and 

Zt{x) = In |t — x| — In |a;|. 

A straightforward computation shows that 

Zt{x) = Sgn(-t) lim / -l[min(x-t,x),max(a;-t,x)](s)rfs, 

SO that, in the space of distributions we get 

Zt = -PV(1/-) * h 
where PV denotes the Cauchy principal value. Thus 

m) = -pv(i7-)(0^(0 

= -(-z7rsgn(0)(-^(e-^«*-l)) 

With f{z) = —-liZt{—z) — l2lt{—z), we obtain 

e*«* - 1 

As in (a) we conclude that fr —>■ f in L". proposition 12.11 implies that F is (7 + l/a)- 
localisable at all m G M with local form Y^ = -liZa + hLa, since M is symmetric. ■ 

We give examples to illustrate propositions 13.11 and 13.21 



f{O = ^7^iih>0 + lh<0) 



Example 3.3 Let | < a < 2 and let M be an a-stable symmetric random measure on M 
with control measure C. Let 

( (x < 0) 

g(x) = I X^l^ (0 < X < 1) . 

\ x-^l^ [x > 1) 
The stationary process defined by 

Y{t) = f g{t- x)M{dx) {t e M) 

is (1/6 + l/a)-strongly localisable at all u &M. with local form Y^ = -^o,i/6+i/a,i,o- 

Proof. We apply proposition 13.11 case (a) with a G (|,2]. The function g satisfies the 
assumptions with •y = ^, Cq = 1, Cq = and a = 1. m 

To verify condition fl3.6p of proposition I3.2[ one needs to check that g G L"(M) and 
also that ^"^^^^^(0 ~ ^^'^ is the Fourier transform of a function in L^(]R) for some a, 7,p 
in the admissible ranges. For this purpose, one may for instance apply classical theorems 
such as in P^ Theorems 82-84]. We give below an example that uses a direct approach. 

Example 3.4 For 1 < a < 2 let M be an a-stable symmetric random measure on M with 
control measure C. Let g be defined by its Fourier transform 



diO - -^ lcl-7-1 



m<K 



where 7 G (— ^,| — ^) ^ (~1)0)- Then g G L"(R) and the moving average process 

Y(t) = [git- x)M{dx) [t G M) 



is well-defined and a-localisable at all u & M., with local form Y^ = La^^+i/a,b,b, where 
6 = -l/(2r(7 + 1) cos(7r(7 + 1/2))). 

Proof. Taking K{^) = |.^|^/^l[_i,i](^) with / = —1 and a = | in (13.61) gives g. To check that 
K G LP{R) for all p > 1, note that K is continuous (in fact C°°) and that \K{x)\ < C|x|"^ 
for all X. Then Y{t) will be well-defined if g is in L"(]R). To verify this, one computes the 
inverse Fourier transform of ^, to get g{x) = 2(7 + l)|x|'^ T , |f |~'^~^ sinvdv — 2x^^ sinx. 
By proposition 13.2( a). Y is a-localisable at all n G M with the local form as stated. ■ 

The approach of this example may be be used for general classes of functions g. 

4 Multistable moving average processes 

In [S], localisability is used to define multistable processes, that is processes which at 
each point t G M have an a(t)-stable random process as their local form, where a{t) is a 
sufficiently smooth function ranging in (0,2). Thus such processes "look locally like" a 
stable process at each t but with differing stability indices as time evolves. 

10 



Before we recall how this was done in [5], we note briefly that "stable-like" processes 
have been defined and studied in [8]. These stable-like processes are Markov jump pro- 
cesses, and are, in a sense, "localisable" , but with localisability defined by the requirement 
that they are solutions of an order a{x) fractional stochastic differential equations. See 
Theorem 2.1 in [8], which shows that the local form of sample paths is considered rather 
than of the limiting process. Another essential difference is that stable-like processes are 
Markov, whereas, in general, multistable ones, as defined below, are not. In fact, formula 
(14.61) . where a Poisson process element Y is independent of t but is raised to a power that 
involves t means that our processes are "far" from Markov. 

We now come back to our multistable processes. One route to defining such processes 
is to rewrite stable integrals as countable sums over Poisson processes. We recall briefly 
how this can be done, see [5] for fuller details. Let {E,S,m) be a a-finite measure space 
and let 11 be a Poisson process on E x M. with mean measure m x C. Thus 11 is a 
random countable subset of £" x M such that, writing N{A) for the number of points in 
a measurable A G E x M., the random variable N{A) has a Poisson distribution of mean 
(m X jC){A) with N{Ai), . . . , N{An) independent for disjoint Ai,. . . ,An C E x R, see 
[6]. In the case of constant a, with M a symmetric a-stable random measure on E with 
control measure m, one has, for j ^ Ta ([101 Section 3.12]), 

f f{x)M{dx) = c{a) J2 /(X)Y<-i/"> (0<a<2), (4.1) 

^ (x,Y)en 

where 

c(a) = {2a-^T{l - a) cosi^na))'^^'^ , (4.2) 

and a^^^ = sign(a)|a|*. 

Now define the random field 

X{t,v)= J2 f{t,v,X)Y<-'/'^(^>. (4.3) 

(x,Y)en 

Under certain conditions the "diagonal" process X{t,t) gives rise to a multistable process 
with varying a of the form 

Y{t)^X{t,t)= Y. /(t,t,X)Y<-i/"W>. (4.4) 

(X,Y)Gn 

Theorem 5.2 of [S] gives conditions on / that ensure that Y is localisable (or strongly 
localisable) with Y^ = X'^{-,u) at a given u, provided X{-,u) is itself localisable (resp. 
strongly localisable) at u. These conditions simplify very considerably in the moving 
average case, taking E = R and m = C with f{t,v,x) = g{x — t). Our next theorem 
restates [5[ Theorem 5.2] in this specific situation. 

We need first to define a quasinorm on certain spaces of measurable functions on E. 
For < a < 6 < 2 let 

^a,b = J^a,b{E,S,'m) = {/ : / is m-measurable with ||/||a,6 < oo} 

where 



a,b 



\f{x)\''m{dx) ) + ( / \f{x)\^m{dx) ) . (4.5) 
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Theorem 4.1 (Multistable moving average processes) Let U be a closed interval with u 
an interior point. Let a : U —* {a,b) G (0, 2) satisfy 

\a{v) — a{u)\ < ki\v — u\''^ {v G U) 

where < t] <1. Let g G J-'a,b, o-nd define 

Y{t)= Y^ ^(X-t)Y<-i/°W> (teM). (4.6) 

(X,Y)Gn 



Assume that g satisfies 

g{r{t + z)) - g{rz) 



lim 



rT 



-h{t,z] 



a(u) 

dz = (4.7) 



for jointly measurable functions h{t, ■) G J^a{u), where < 7 + l/a{u) < f] < 1. Then Y 
is (7 + l/a{u))-localisable at u with local form Y^ = {J h{t,z)Ma{u){dz) : t G M}, where 
Ma[u) is the symmetric a{u) -stable measure with control measure C and skewness 0. 
Suppose further that 7 > and for h sufficiently small 

\\g{h-x)-gi-x)\U<c\h\^-^'/"^^\ 

Then Y has a continuous version and is strongly (7 + 1 / a{u)) -localisable at u with local 
form Y^ = {J h{t, z)Ma{u){dz) : t G M} under either of the following additional conditions: 

(i) < a{u) < 1 and g is bounded 

{ii) 1 < a{u) < 2 and a is continuously differentiable on U with 

\a' iy) — a' {w)\ < ki\v — wY' {v,weU), 
where l/a{u) < rj < 1. 

Proof. Taking 

X{t,v)= J2 ^(X-t)Y<-i/"(^> (t,wGM). (4.8) 

(x,Y)gn 

this theorem is essentially a restatement of O Theorem 5.2] in the special case of i? = R 
and m = C with f(t,v,x) = g{x — t) in (14. 3p . Since f(t,v,x) no longer depends on v 
most of the conditions in [5i Theorem 5.2] are trivially satisfied and we conclude that 
FJ = X'^{-,u), noting that X{-,u) is (7 + l/a;(M))-localisable (or strongly localisable) with 
the local form given by propositions 12.11 or 12.41 ■ 

It is curious that neither cases (i) or (ii) address localisability if a{u) = 1. This goes 
back to the proof of [3, Theorem 5.2] where different approaches are used in the two cases. 
For a{u) < 1 the proof uses that the sum (4.8) is absolutely convergent almost surely, 
whereas for a{u) > 1 we need to find a number p such that < p < a{t) for t near u with 
r^p > 1 to enable us to apply Kolmogorov's criterion to certain increments. 

Corollary 4.2 Let U,a and g be as in Theorem \4.1\ Then the same conclusion holds if 
Y{t) m ^^ IS replaced by Y{t) = a{t) E(x,y)en di^ - t)Y<-i/"(*)> (t G M), where a is 
a non-zero function of Holder exponent rj > h. 
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Proof. This follows easily in just the same way as proposition 2.2 of [5]. ■ 

We may apply this theorem to get a multistable version of the reverse Ornstein- 
Uhlenbeck process considered in Section [21 

Proposition 4.3 (Multistable reverse Ornstein- Uhlenbeck process) Let A > and a : 

M -^ (1, 2) be continuously dijjerentiable. Let 

Y{t)= Yl exp(-A(X-t))Y<-^/"W> (tER). 
(x,Y)en,x>t 

Then Y is 1 / a{u) -localisable at all u eM with Y^ = c{a{u))^^La{u), where La is a-stable 
Levy motion. 

Proof. Taking g{x) = l[o,oo)(2^) exp(— Aa;) and h(t, z) = — 1[_( o](^) for t > (and a similar 
formula for t < 0) with 7 = 0, localisability follows from Theorem 14 . 1 1 with the limit ( 14. 7p 
being checked just as in proposition I2.2[ b 

Theorem 14 . 1 1 applies in particular to functions g satisfying the conditions of proposition 
13. 1[ Thus, for instance, the moving averages of Examples 13.31 and 13.41 admit multistable 
versions. The process of Example 13.31 is strongly 7 + l/a{u) localisable at u whenever a 
verifies condition (ii). 

5 Path synthesis and numerical experiments 

We address here the issue of path simulation. In the previous sections, we have considered 
two kinds of stochastic processes: moving average stable ones, that are stationary, and 
their multistable versions, which typically are not, nor have stationary increments. Our 
simulation method for the moving average stable processes is based on that presented in 
|12j . There, the authors propose an efficient algorithm for synthesizing paths of linear 
fractional stable motion. In fact, this algorithm really builds traces of the increments of 
linear fractional stable motion. These increments form a stationary process, an essential 
feature for the algorithm to work. It is straightforward to modify it to synthesize any 
stationary stable process which possesses an integral representation. In addition, we are 
able to obtain bounds on the approximation error measured in the a-norm, and thus on 
the r*'' moments for r < a, as shown below. 

For non (increment) stationary processes, like multistable processes, a possibility 
would be to use the general method proposed in [1]. It allows to synthesize (fractional) 
fields defined by integration of a deterministic kernel with respect to a random infinitely 
divisible measure. When the control measure is finite, the idea is to approximate the 
integral with a generalized shot noise series. In this situation, a bound on the L^ norm 
of the error is obtained for appropriate r. In the case of infinite control measure, one 
needs to deal with the points "far from the origin" through a normal approximation. 
This second approximation maybe controlled through Berry-Esseen bounds which lead to 
a convergence in law. Thus the overall error when the control measure is infinite may 
only be assessed in law, and not in the stronger L^ norm. 

Although the method of [IJ may be used for the synthesis of multistable processes, 
we will rather take advantage here of the particular structure of our processes: being 
localisable, they are by definition tangent, at each point, to an increment stationary 
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process. Thus we may simulate them by "gluing" together in a appropriate way paths of 
their tangent processes, which are themselves synthesized through the simpler procedure 
of [12] . This approach allows in addition to control the error in the U norm, rather than 
only in law (since we are in an infinite control measure case). 

We briefly present in the next subsection the main ingredients of the method. We 
then give bounds estimating the errors entailed by the numeric approximation, in the 
case where the process is localisable. Finally, we display graphs of localisable moving 
average processes obtained with this synthesis scheme. 

5.1 Simulation of stable moving averages 

Let Y = {Y{t),t G M} be the process defined by (12. 3p . To synthesize a path Y{k),k = 
1, ...,N,N G N, oiY, the usual (Euler) method consists in approximating the integral by 
a Riemann sum. Two parameters tune the precision of the method: the discretization 
step uj and the cut-off value for the integral fi. The idea in [12] is to use the fast Fourier 
transform for an efficient computation of the Riemann sum. More precisely, let 

Y{k) = f g{k- s)dM{s) = - f g{s)dM{k - s). 

Jr Jr 

Let cu, r2 G N and 

. ^ ujQ . 

YMk)= T 9i^-^)ZaA^k-j) + y2g{^)Z^,^{uk-j), (5.1) 

j=-uin+i j=i 

where Za^j) = M(^^) — M(-) are i.i.d. a-stable symmetric random variables. Let 
Za{j) denote a sequence of normalised i.i.d a-stable symmetric random variables. Then 
one has the equality in law: {Za^uiij),j G Z} = {uj^^^°'Za{j),j G Z}. One may thus write: 

2a;n 

Y^Ak) = J2(^^ij)Za{oo{k + Q) -j), 
j=i 

where 

i^-'^"9i^-^) forjG{l,...,cufi} 

a;-i/°t/(i - Q) for j G {uQ + 1, ..., 2uQ}. 

For n G Z, let 

2ujn 

i=i 

Then {Y^,nik),k = 1,...,A^} has the same law as {W{u{k + n)),k = 1,...,7V}. But W 
is the convolution product of the sequences a^ and Za- As such, it may be be efficiently 
computed through a fast Fourier transform. See [12] for more details. 



•^^0) - ^ , ,-l/a^(l _ 



5.2 Estimation of the approximation error 

When the moving average process is localisable, or more precisely when the conditions of 
proposition 13.11 are satisfied, it is easy to assess the performances of the above synthesis 
method. 
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The following proposition gives a bound on the approximation error in the a— norm. 
Recall that the a— norm (defined in (12. ip ) is just the scale factor of the random variable, 
and is thus independent of the integral representation that is used. In addition, it us, up 
to a constant depending on r and a, an upper bound on moments of order < r < a. 

Proposition 5.1 Let Y be defined by Ii2.3\) . and let Y^^^q be its approximation defined in 
( 15. ij) . Assume g satisfies the conditions of proposition [X7l Then, for all u!,Q E N and 
k E "L with to > -, one has 

Err:=\\Y{k)-Y^,n{k)h<Al{^ 
where 

A^n = - r^^V^ — ^+ / \9(s)rds+ \g(sTds 

Proof. By stationarity and independence of the increments of Levy motion, one gets: 
Err= J2 \9{- ) - 9{s)rds + J2 \9{^) - gisTds 

/—SI /•+00 

\g{sTds+ \gis)rds. (5.2) 

By assumption, for almost all s G M, \g{s + h) — g{s)\ < c\h\"'\s\'^^'^ when < h < rj. 
Recall that u > -. A change of variables yields 

rJ/'^ 7 — 1 /'^/'^ 7 — 1 
/ \g{^- )-g{sTds< c'^lsn^- |(^~")"ds 

and thus 

Err < V / c'^lsri- |(^-'^)"rfs + V / c"|sr"|-|(^-")"ds 

/ — SI r+oo 

\gis)rds+ / IgisTds. 
-oo Jfl 

Rearranging terms: 



Err < 



^ 


c" 1 




(l + aa)cji+'y" 




.-Q 


+ 


/ \9{sTds- 




J —CO 


= 


A^n- 



\j=-ujn+i j=i / 



which is the stated result. ■ 
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Corollary 5.2 Under the conditions of proposition l5A\ \\Y{k) — Y^,^{k)\\a -^ when 
{uj, Q) tends to infinity. 

If in addition g[x) < C\x\^^ when |x| — >■ oo for some C > and P > -, then: 

\\Yik) - Y^Mkm < K (cu-i--^ + fii-"^) (5.3) 

where K is a constant independent of k,uj,Q. 

Proof. Since g satisfies the assumptions of proposition 13. ![ a > 7 + ^. As a consequence, 
the sum in the first term of A^^q converges when {u, Q) tends to infinity. The first 
statement then follows from the facts that 07 + 1 > and g G J^a- The second part 
follows by making the obvious estimates. ■ 

The significance of (15. 3p is that it allows us to tune uj and Q to obtain an optimal 
approximation, provided a bound on the decay of g at infinity is known: optimal pairs 
{uj, Q) are those for which the two terms in (15. 3p are of the same order of magnitude. More 
precisely, if the value of (3 is sharp, the order of decay of the error will be maximal when 

— l — a.'y -. 

Q = UJ i-"/3 . Note that the exponent "J^*? is always positive, as expected. Intuitively, uj 
is related to the regularity of g (irregular g requires larger uj), while Q is linked with the 
rate of decay of g at infinity. 

For concreteness, let us apply these results to some specific processes: 

Example 5.3 (reverse Ornstein-Uhlenbeck process) LetY be the reverse Ornstein- Uhlenbeck 
process defined in proposition \2.^ When a > 1, we may apply proposition I5.il with 
g{x) = exp(x)l(a; < 0),'y = 0, a = 1, c = 2, rj = 1. One gets, for uj > 1, il > 1, 

^a+l ( ^^ 




However, we may obtain a more precise bound on the approximation error, valid for 
any a G (0, 2), by using (15.21) directly: 



„— aSl \ 1 p— qSI 



r (A:) - Y^^^{k)\\l < —— — ^— ^^ + 

1 + a \ e°'i^ ~ 1 / ^ « 

When {uj,VL) -^ +00, Err < 0{-^) + (9(e~"^), which is better than A^,n above when 
a > 1. 

We note finally that the optimal choice for {uj, VL) is here Vt = ln{uj), which is consistent 
with the fact that the f3 in Corollay 15.21 may be chosen arbitrarily large. 

Example 5.4 (linear fractional stable noise) Let < a < 2 and let M be an a-stable 
symmetric random measure on M with control measure C. Let: 

g{x) = (x)f ~^/" -{x- 1)+"^/" 

and 

Y{t) = f g{t- x)M{dx) {t G M) 



Applying the analysis above with 7 = if — -, a = l,c = 2|7|,r7 = 1 one gets, for u; > 1, 

n>i, 

Oa+l \ II ]_\a '^^ -I -I r+00 
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When {uj, il) -^ +cxd, 

1 /"+00 

^"^ Jn 



C(— -) + c)(^i+°(^-i/"-i)^ 



This process is the one considered in [12]. Here we reach a conclusion similar to 
Theorem 2.1], which yields the same order of magnitude for the error when [u, Q) — > +oo. 
Extensive tests are conducted in [12] to choose the best values for {uj,Q). The criterion 
for optimizing these parameters is to test how an estimation method for H performs on 
synthesized traces. Here we adopt a different approach based on Corollary 15.21 optimal 
pairs {u, Q) are those for which (15.31) is minimized. Since the value oij3 = l — H + l/a 
is sharp here, one gets Q = uj^-h . It is interesting to note that the exponent H/{1 — H) 
depends only on the scaling factor H and not on a, and that it may be larger or smaller 
than one depending on the value of H. We do not have an explanation for this fact nor 
for the reason why H = 1/2 plays a special role. 

Example 5.5 As a final illustration, we consider the process of Example \3.3[ . Withj = ^, 
a = 1,0= l,rj = 1, one gets, for lj > 1, Q > 1, 

ct+6 

Again, the value of /5 = 5/6 is sharp, and the optimal choice is to set i7 = cjsa-s. 
Since {a + 6)/{5a — 6) > 1, fi is larger than a; in this case, in contrast to the reverse 
Ornstein-Uhlenbeck process: it is the decay at infinity of the kernel that dictates the 
parameters here, while it was the regularity that mattered in the case of the reverse 
Ornstein-Uhlenbeck process. 

5.3 Numerical experiments 

We display in figure [T] traces of: 

• moving average stable processes: the reverse Ornstein-Uhlenbeck process (figures 
l(e),(f)), and the processes of Examples 13.31 (figure 1(c)) and 13.41 (figure 1(a)). In 
each case, a = 1.8. Some of the relevant features of the processes of Examples 
13.31 and 13.41 seem to appear more clearly when one integrates them. Thus integral 
versions are displayed in the right-hand part of the corresponding graphs, figures 
l(b),(d). 

• a multistable version of the reverse Ornstein-Uhlenbeck process, using the theory 
developed in Section H] (figures 1(g), (h)). Since these processes are localisable, one 
may obtain paths by computing first stable versions with all values assumed by a, 
and then "gluing" these tangent processes together as appropriate. More precisely, 
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assume we want to obtain the values of a multistable process S at the discrete 
points (ti,...,t„). At each ti, S is tangent to a stable process denoted Si. We 
first synthesize the n stable processes Si with the method just described, all with 
the same random seed. The multistable process S is then obtained by setting 
S(ti) = Si(ti),i = 1, . . . ,n. Two graphs are displayed for the multistable process: 
in figure 1(g) the graphs are as explained above. In figure 1(h) each "line" of the 
random field {i.e. the process obtained for a fixed value of a) is renormalized so 
that it ranges between -1 and 1, prior to building the multistable process by gluing 
the paths as appropriate. This renormalization may be justified using Corollary 14. 2[ 

The parameters are as follows: 

• Process of Example 13.41 uj = 5000, Q = 877, N = 2000. The approximation error 
Err is bounded by 2.172. 

• Process of Example E^ u = 104, Q = 175504, A^ = 2000. The term AI[^ is equal 
to 0.074. 

• Reverse Ornstein-Uhlenbeck process with \ = 1 : u = 512, Q = 7,N = 7392. The 
term AJ'^ is equal to 0.0018. 

• Reverse Ornstein-Uhlenbeck process with A = 0.01 : u = 256, fi = 800, A^ = 7392. 
The term AJ^ is equal to 0.0032. 

• Multistable reverse Ornstein-Uhlenbeck process: A = 0.01, a; = 256,^2 = 800, A^ = 
7392. The a function is the logistic function starting from 1.2 and ending at 1.85. 
More precisely, we take: a(t) = 1.2 + , ^f^,^ ^r,^^^ ^ where A^ is the number of 

points and t ranges from 1 to A^ (the graph of a{t) is plotted in figure 1(h). Thus, 
one expects to see large jumps at the beginning of the paths and smaller ones at 
the end. Note that we do not have any results concerning the approximation error 
for these non-stationary processes. 

The value of Q in all cases is adjusted so that the pair {uj,Q) is approximately "optimal" 
as described in the preceding subsection (optimality is not guaranteed for the multistable 
processes. Nevertheless, since the relation between u and Q does not depend on a for the 
reverse Ornstein-Uhlenbeck process, it holds in this case). 

The function g of example 13.41 cannot be treated using corollary 15.21 nor proposition 
15.11 since g does not satisfy the conditions of proposition 13. 1[ However, it is possible 
to estimate Err directly. Since \g{x + h) — g{x)\ < 2|/i|(|x|~2l|^|^]^ + |x|~^l|^l>i) and 
\g{x)\ < ^^, one gets: 

)a+2 ^^ 



^ „ 2"+2 a 1 8 1 

Err'' < -^-^ + 



l + aa-lu;^~2 a - 1 fi"- 



2-a 



The asymptotic optimal relation between u and Q is thus ^2 = ^^^{q-i) = tu^'^^^. The 
values in our simulation are slightly different since they are chosen to optimize the actual 
expression with a finite a;. 
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Finally, we stress that the same random seed {i.e. the same underlying stable M{dx)) 
has been used for all simulations, for easy comparison. Thus, for instance, the jumps 
appear at precisely the same locations in each graph. Notice in particular the ranges 
assumed by the different processes. 

The differences between the graphs of the processes of Examples 13.31 13.41 and the 
reverse Ornstein-Uhlenbeck process are easily interpreted by examining the three kernels: 
the kernel of the process of Example 13.41 diverges at 0, thus putting more emphasis on 
strong jumps, as seen on the picture, with more jaggy curves and an "antipersistent" 
behaviour. 



The kernel of the process of Example 13.31 in contrast, is smooth at the 
origin. In addition, it has a slow decay. These features result in an overall smoother 
appearance and allow "trends" to appear in the paths. Finally, the kernel of the reverse 
Ornstein-Uhlenbeck process has a decay controlled by A. For "large" A (here, A = 1), little 
averaging is done, and the resulting path is very irregular. For "small" A (here, A = 0.01), 
the kernel decays slowly and the paths look smoother (recall that, in the Gaussian case, 
the Ornstein-Uhlenbeck tends in distribution to white noise when A tends to infinity, and 
to Brownian motion when A tends to 0). 
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Figure 1: Paths of locahsable processes, (a) The process in Example 13.41 and (b) the 
integrated version, (c) The process in Example 13.31 and (d) the integrated version, (e) 
Reverse Ornstein-Uhlenbeck processes with A = 1, and (f) A = 0.01. (g) A multistable 
reverse Ornstein-Uhlenbeck process with A = 0.01 and (h) the renormalised version along 
with a{t). 
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